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We simulate the ordering of vortices and its effects on the critical current in superconductors with 
varied vortex-vortex interaction strength and varied pinning strengths for a two-dimensional system. 
For strong pinning the vortex lattice is always disordered and the critical depinning force only weakly 
increases with decreasing vortex- vortex interactions. For weak pinning the vortex lattice is defect 
free until the vortex-vortex interactions have been reduced to a low value, when defects begin to 
appear with a simultaneous rapid increase in the critical depinning force. In each case the depinning 
force shows a maximum for non-interacting vortices. The relative height of the peak increases and 
the peak width decreases for decreasing pinning strength in excellent agreement with experimental 
trends associated with the peak effect. We show that scaling relations exist between the distance 
between defects in the vortex lattice and the critical depinning force. 



I. INTRODUCTION 

A large class of condensed matter systems can be rep- 
resented as an elastic lattice interacting with a quenched 
random substrate. Vortices in type II superconductors, 
which can be pinned by random defect sites in the ma- 
terial, are a particularly ideal example of such a system 
since the substrate and lattice interactions can be read- 
ily tuned experimentally. In this case, a rich variety of 
static and dynamical phases can occur due 

to the competition between the repulsive vortex-vortex 
interactions, which tend to order the vortex lattice, and 
the attractive vortex-pin interactions, which tend to dis- 
order the vortex lattice. Unlike the early studies of col- 
lective pinning of a defect-free elastic lattice, recent ex- 
perimental and theoretical work point to the importance 
of topological defects in pinned vortex matter. If the 
pinning in the sample is sufficiently strong, the vortex 
lattice will be highly defected; conversely, in a sample 
with weak pinning the vortex lattice can be relatively or- 
dered and contain only a few defects. In contrast with the 
purely elastic model, however, the relationship between 
the critical force and the density of flux line defects is not 
known. Moreover, applications of superconducting ma- 
terials require that the critical current be high, and so a 
key question is how the number of defects in the vortex 
lattice can affect the critical current. 

In a given sample, i.e., with a given realization of 
quenched disorder, the relative strength of the pinning 
can be increased by lowering the vortex-vortex interac- 
tion strength, i.e., by softening the lattice. Vortex- vortex 
interactions become weak at low applied fields, when vor- 
tices are dilute, as well as very near H C 2 and T c , where 
the induction of the vortices overlap significantly. In the 
latter regime one encounters the peak effect: a rapid in- 
crease in the critical current before it collapses to zero at 



the superconducting - normal phase boundary. At these 
two extreme field values, the pinning interaction domi- 
nates and the vortex lattice is expected to be highly dis- 
ordered. Indeed, direct evidence for vortex lattice disor- 
dering through the peak effect regime has been observed 
in neutron scattering measurements ^| . A reentrant dis- 
ordering has been found in the low field regime where the 
vortex lattice again softens Q . Bitter decoration experi- 
ments have also provided evidence for the disordering of 
the vortex lattice at low fields Jl0| . In each of these cases 
the critical current J c is high due to the dominance of 
pinning. 

Experiments on superconducting samples with differ- 
ent amounts of pinning have revealed some systematic 
trends. The peak effect is much sharper and more pro- 
nounced in cleaner samples [pp|Jll]| , while in samples with 
stronger pinning, the width of the peak region increases 
and the relative valley-to-peak height decreases, although 
the overall critical current is higher In samples 

with the strongest pinning, J c is high through most of 
the applied field range, and the peak is completely ab- 
sent. A possible interpretation of this behavior connects 
vortex lattice defects with high critical currents. A de- 
tailed understanding of the observed effects as a possible 
transition/crossover between regimes with low and high 
critical currents, i.e., between, say, a Bragg glass- like or- 
dered phase with few defects and a Vortex glass/ pinned 
liquid-like disordered vortex phase is still lacking. A key 
reason for this uncertainty is the absence of a systematic 
connection between the number of defects in the vortex 
lattice, the softness of the lattice, and the critical current. 

In this work, we directly examine the effect of vor- 
tex lattice softness on the critical current through a se- 
ries of numerical simulations of samples with different 
pinning strengths. For strong pinning the vortex lat- 
tice is highly defective for the entire range of vortex- 
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vortex interactions investigated. The critical depinning 
force, which is proportional to the critical current, only 
weakly increases with decreasing vortex-vortex interac- 
tion strength, reaching a maximum for non-interacting 
vortices. For the weakest pinning strengths, on the other 
hand, the lattice is almost defect-free over a large range of 
vortex- vortex interaction strengths. When defects begin 
to appear in the softest vortex lattices, the critical depin- 
ning force increases sharply from its low value to a peak 
at zero vortex-vortex interaction. The behavior for inter- 
mediate pinning is intermediate between these extremes. 
We find that the relative height of the critical current 
peak increases for decreasing pinning while the width of 
the peak decreases. Thus, the most pronounced peaks 
occur for the weakest pinning. These results are in excel- 
lent agreement with the experimental trends described 
above. We discuss how these results can be connected to 
the behavior of the peak effect in superconductors with 
varying pinning strength. 

II. SIMULATION 

We consider a 2D slice of a system of superconducting 
vortices interacting with a random pinning background. 
The applied magnetic field H = Hi is perpendicular to 
our sample, and we use periodic boundary conditions in 
x and y. The T = overdamped equation of motion for 
a vortex i is 

where Vj is the velocity of vortex i and r\ = 1 is the 
damping coefficient. The total force on vortex i from 
the other vortices is /™ = A v fo K i(\ r i ~ r j|A) r yj 

where r< is the position of vortex i, A is the penetration 
depth, fo — <&g/87r 2 A 3 , the prefactor A v is used to vary 
the vortex lattice softness, and K\ is the modified Bessel 
function. The pinning f^ p is modeled as randomly placed 
attractive parabolic wells of radius r p — 0.15A with /" p — 
(f p /r p )(\n - rjWj) Q(r p - |r, - r^|)f^, where r[ p) is 
the location of pin k, f p is the maximum pinning force, 
which is varied from 0.1/o to 3.0/o, 6 is the Heavside 
step function and = (r* — rj. p ')/|ri — |. The pin 
density is n p = 3.0/ A 2 and the vortex density is n v = 
0.75/A 2 . We simulate a 36A x 36A system containing 
N v = 864 vortices and N p — 3887 pins. We initialize 
the vortex positions by performing simulated annealing, 
starting from a high temperature and slowly cooling to 
T = 0. This method of preparing the lattice is similar to 
field cooled experiments. To identify the depinning force 
f c we apply a slowly increasing uniform driving force fd 
on the vortices in the x direction, which would correspond 
to a Lorentz force from an applied current J = Jy. For 
each drive increment we measure the average vortex 




(e) (f) 
FIG. 1. The Delaunay triangulation of the vortex lattice 
for decreasing vortex-vortex interaction strength A v in a sam- 
ple with f p — 0.25/o • Dark circles indicate 5-fold coordinated 
vortices; open circles indicate 7-fold coordinated vortices, (a) 
Av = 4.0, (b) 3.0, (c) 2.0, (d) 1.0, (e) 0.75 and (f) 0.50. The 
vortex lattice is relatively ordered in (a). In (b) more defects 
appear, (c) and (d) show different domains of vortex orienta- 
tion. In (e) and (f) the vortex lattice is highly disordered. 

velocity in the direction of drive, V x — (1/N V ) v x- 
The fd versus V x curve corresponds experimentally to a 
V(I) curve. The depinning force f c is defined as the drive 
at which K. > 0.03. 



III. VORTEX ORDER AND PINNING FOR 
VARIED VORTEX- VORTEX INTERACTION 
STRENGTH 

We first consider the effect of the vortex lattice softness 
on the stationary vortex lattice. In Fig. [j] we show the 
Delaunay triangulation for a system with f p = 0.25/o for 
decreasing vortex- vortex interaction A v , after the 
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FIG. 2. (a) The critical depinning force f c versus A v for a 
sample with f p = 0.25/o- (b) The fraction of six-fold coordi- 
nated vortices P@ versus A v . As the vortex- vortex interaction 
is lowered f c increases while Pq decreases. 

lattice has been annealed and with no driving force ap- 
plied. Defect sites in the vortex lattice are indicated 
by circles. In Fig. 0(a) for A v = 4.0, the vortex lat- 
tice contains no defects. In Fig. 0(b) for A v — 3.0 a 
small number of 5-7 defect pairs appear. In Fig. 0(c), 
at A v = 2.0 the vortex lattice is considerably disordered 
but some domains of order are still present. For softer 
vortex lattices in Fig. 0(d-f), the vortex lattice becomes 
progressively more disordered by the underlying random 
pinning. Since the simulation is for a two-dimensional 
system, we expect the system to be defective for any 
pinning strength. The defect-free situation likely results 
from the finite size of the sample. 

By applying a transport current to the annealed lat- 
tices, we determine the critical depinning force f c that 
must be applied before the vortices begin to move. In 
Fig. ||(a) we plot f c versus A v for a sample with f p — 
0.25/o. For comparison, we measure the amount of order 
in the lattice before depinning using the Delaunay tri- 
angulation, and in Fig. ||(b) we plot Pq, the fraction of 
six-fold coordinated vortices, versus A v . Here Pq = 1.0 
indicates a perfect triangular lattice. As A v is lowered f c 
increases and simultaneously the order in the lattice, Pq, 
decreases. The maximum value of f c occurs at A v ~ 
which coincides with the minimum value of Pq. Thus the 
softer lattices with low values of A v and large amounts 
of disorder are more strongly pinned than stiffcr, more 
ordered lattices. 




FIG. 3. Scaled critical depinning force and Pq versus 
A v for samples with different pinning strengths. Open cir- 
cles: f p = 0.10/o; filled squares: f p = 0.25/o; open dia- 
monds: f p — 0.75/o; filled triangles: f p = 1.5/o; plus signs: 
fp — 3.0/o. (a) The plot of f c /f P versus A v shows that 
the peak at low A v becomes sharper for weakening pinning 
strength f v . (b) The corresponding Pa values show that for 
large f p = 3.0 the vortex lattice is disordered over the entire 
range of A v . 



IV. HEIGHT AND WIDTH OF CRITICAL 
DEPINNING FORCE PEAK 

The shape and magnitude of the critical current peak 
at A v — are affected by the strength of the pinning in 
the sample. To demonstrate this, in Fig. ||(a) we plot 
fc/ ' fp versus A v . Scaling the curves with f p in this way 
causes all of the curves to approach f c /f p = 1 at A v = 0. 
We find that the relative height of the peak increases 
as the pinning becomes weaker. For example, compar- 
ing f p — 0.25 (bottom curve) with f p = 3.0 (top curve), 
the critical current f c increases by a factor of 8 from 
f c = 0.1 at A v = 6.0 to f c = 0.8 at A v = 0.01 in the 
weakly pinned sample, whereas in the strongly pinned 
sample, f c increases only by a factor of 2 from f c = 0.38 
at A v = 6.0 to f c = 0.8 at A v = 0.01. The peak not 
only becomes higher for weaker pinning, but also it be- 
comes much narrower as can be seen by comparing the 
widths of the peaks in Fig. ||(a). As shown in Fig. ||(b), 
in the strongly pinned samples such as f p = 3.0/o with 
broad low peaks, P§ is low over the entire range of A v , 
indicating that the vortex lattice is always strongly dis- 
ordered by the pinning. In contrast, for weakly pinned 
samples such as f p = 0.25/o, P% drops a large amount 
from Pq w 1 at large A v , indicating a nearly perfect 
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FIG. 4. Squares: the relative 

height H = f c (A v = 0.01) // C (A„ = 6.0) of the peak in f c 
for samples with different pinning strength f p . Circles: The 
width of the peak dW at half-maximum as a function of f p . 
The peak becomes sharper and more pronounced for weaker 
pinning. 

lattice, to Pq « 0.4 at the lowest A v , indicating a large 
amount of disorder in the lattice. 

In Fig. [| we show explicitly how the relative peak 
height, H = f c (A v = 0.01)// C (A„ = 6.0), and the width 
at half maximum, dW, vary with pinning strength f p . 
Here we see that the relative strength and sharpness of 
the peak increase with weaker pinning. The width dW 
depends linearly on f p , while the relative height H shows 
a nonlinear increase that can not be fitted by a simple 
functional form. 

The results in Fig. || and Fig. |] suggest that in clean 
superconductors, where the vortex lattice softness has 
a much more significant impact on the effectiveness of 
the pinning, the critical current enhancements should be 
both sharper and of larger relative height than in strongly 
pinned superconductors. This behavior agrees with the 
results of recent experiments |||| on the peak effect. Our 
results also suggest that the peak can be made arbitrar- 
ily sharp simply by adjusting the disorder strength, and 
that such large peaks in the critical current must be ac- 
companied by the proliferation of defects. 

V. CRITICAL DEPINNING FORCE 
DEPENDENCE ON DEFECT DENSITY 

Using the data obtained by varying A v , we can de- 
termine the relationship between the critical current and 
the density of defects. In Fig. |](a) we plot f c /fp versus 
Pq, showing that the curves for samples with different 




FIG. 5. (a) Scaled critical depinning force f c / f p versus 
Pa. (b) fc/fp versus the distance between defects d<j showing 
a collapse for dd < 3.0 coinciding with a rapid increase in f c . 

pinning strengths roughly collapse. At lower values of 
Pq when the vortex lattice becomes defected, the critical 
current increases rapidly with the onset of defects. In 
Fig. |(b) we plot f c /f p versus the distance between vor- 
tex lattice defects, dd = 1/y/TH — Pq. In samples with 
weak pinning there is a range of high A v values that pro- 
duce well ordered vortex lattices. In this regime, when 
there are very few defects in the vortex lattice which are 
far apart, we find that the f c /f P values do not collapse, 
and that the rate at which f c /fp changes with dd is slow. 
For dd < 3.0, or when the lattice contains approximately 
10% or more defects, all the curves collapse and f c /f P 
increases rapidly with decreasing dd- The behavior of f c 
versus dd suggests that there is a critical distance between 
defects below which J c begins to rapidly increase. 

We can compare the dependence of the critical de- 
pinning force on the defect density with the collective 
pinning results of the Larkin-Ovchinnikov (LO) theory 
p"2[ . In Fig. ||(a) we plot f c as predicted by LO, using 
f c (x fp/A v . Fig. ^|(b) shows f c computed from the dis- 
tance between dislocations, using f c cx f p /dd- The value 
of f c measured in the simulations is shown in Fig. |^(c). 
We compare f c obtained by the LO and dislocation dis- 
tance methods with the actual f c in Fig. ||(c) . The agree- 
ment with the LO prediction is poor, but that with the 
dislocation distance is good out to values A v /f p ~ 5, 
beyond which dd becomes of order the system size. 
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two-dimensional lattices with quenched disorder where 
the lattice interaction can be tuned. 




A/f A/f 



FIG. 6. (a) Shaded symbols: f c /fp as predicted by LO 
theory, (b) Filled symbols: f c /f P computed from the distance 
between dislocations, (c) Open symbols: f c /f P measured in 
the simulations, (d) All three methods of obtaining f c are 
scaled by the value at A v = 0.01, and plotted together. For all 
panels, circles: f p — 0.10/o; squares: f p — 0.25/o; diamonds: 
f p = 0.75/o; triangle up: f p — 1.5/o; triangle left: f p — 3.0/o. 

VI. CONCLUSION 

We have investigated the dependence of the critical 
dcpinning force f c and vortex lattice topological order 
P 6 on the vortex lattice rigidity for different values of 
pinning strength. The study is restricted to the two di- 
mensional case which is appropriate for vortex lattices 
where the longitudinal correlation length is large enough 
to be effectively given by the sample thickness. In all 
cases, f c increases with decreasing vortex-vortex interac- 
tion strength, reaching a maximum for non-interacting 
vortices, as expected. Although for very weak pinning 
and strong vortex-vortex interaction, the lattice is rela- 
tively defect-free, for strong pinning the lattice is defec- 
tive for the entire vortex interaction range investigated. 
These results are in accord with previous studies. Wc 
have, in addition, studied the systematic variation of f c 
with varying pinning and inter-vortex interaction. The 
results show that the pinning force variation departs sig- 
nificantly from the expectations of an elastic picture. The 
increase in f c towards the maximum occurs very rapidly 
with decreasing interaction or increasing pinning. The 
relative height increases while the width decreases with 
decreasing pinning strength, strongly remniscent of the 
peak effect phenomenon. Furthermore, a pseudo-Larkin 
picture where a plastic length, the distance between topo- 
logical defects, replaces the elastic correlation length, 
provides a good account of the variation of f c . Wc also 
find evidence that there is a critical value of the plastic 
length, typically spanning 3-5 vortices on each side, at 
which the rapid crossover occurs. These results should 
be relevant to experiments in very weak pinning (quasi- 
2D) flux lattices, as well as to the general systems of 
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